Eigenvalue repulsion in an effective theory of SU(2) Wilson lines in three dimensions 
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We perform simulations of an effective theory of SU(2) Wilson lines in three dimensions. We 
include a non-perturbative "fuzzy-bag" contribution which is added to the one-loop perturbative 
potential for the Wilson line. We confirm that, at moderately weak coupling, this leads to eigenvalue 
repulsion in a finite region above the deconfining phase transition which shrinks in the extreme weak- 
coupling limit. A non-trivial Z(N) symmetric vacuum arises in the confined phase. 
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I. INTRODUCTION 

In QCD at very high temperature, the pressure is due 
to weakly interacting quasi-particle gluons. (Here, we 
ignore contributions from quarks and focus on the pure 
gauge theory.) It can be calculated from an effective 
theory in three dimensions 1 , 



£ eff = 



^trGl+tr\D t A \ 



i D tr A 



(1) 



see [J, [2| and references therein. Gij is the magnetic field 
strength associated with the spatial components of the 
vector potential A, mo is the Debye mass, and the dots 
represent self-interactions of Aq. The effective theory (JTJ 
is valid as long as fluctuations in Aq are small; that is, 
the expectation value of the Polyakov loop in the original 
four-dimensional theory, which is given by the trace of the 
thermal Wilson line, 

l(x) = ^trL(x) , (2) 

( 7 \ 

L(x) = Z^Vexp ig / dr A (x,t) , (3) 



should be close to one of the N roots of unity, where N is 
the number of colors. We have indicated explicitly that 
in the four-dimensional theory Aq(x,t) depends on Eu- 
clidean time r and that Polyakov loops have to be renor- 
malized to obtain a non-zero continuum limit [3, |1] . The 
renormalization constant Z~ x depends on the represen- 
tation of L, taken here to be the fundamental represen- 
tation; thus, L represents the propagator of an infinitely 
heavy test quark. Also, we always consider normalized 
traces and divide by the dimension of the representation. 
The Polyakov loop represents an order parameter for 
the spontaneous breaking of the global Z(N) center- 
symmetry corresponding to gauge transformations which 



are periodic in r only up to an element of Z(N). In the 
high-temperature deconfined phase, the Polyakov loop 
aquires a non-vanishing expectation value but vanishes 
in the confined phase @. As a consequence of the 
Z(N) symmetry in the original four-dimensional theory, 
when Aq/T is large (of order 1/g), the effective elec- 
tric field in three dimensions is not simply Ei(x) — 
Di(x)A (x) 0,i. 

For two colors, by universality [7J the phase transition 
is of second order § and so {(){T d ) = vanishes continu- 
ously. Hence, in the immediate vicinity of Td at least, the 
Polyakov loop is clearly far from unity. For N = 3, the 
transition is first-order Q and (£} is discontinuous at Td- 
If the SU(3) gluon plasma was perturbative all the way 
down to T+ (in the electric sector, (^)(T+) « 1), then Q 
might have applied even for T just above Td. However, 
lattice measurements [3,11] indicate that {(){T^) ^0.5, 
which is rather far from unity. Moreover, the ratio of 
screening masses defined from two-point correlation func- 
tions of the real and the imaginary part of the Polyakov 
loop, respectively, increases from «3:2 (which is the LO 
result from perturbation theory) at high temperature to 
«3 near Td [101 ]. Also, resummations of perturbation 
theory work very well at high T but appear to fail to 
reproduce the pressure or the entropy density [ll| below 
ss 3Td- Finally, the interaction measure (e — 3p)/T 4 for 
both two [3] and three Q colors is rather large up to 
T w 3T d . 

These observations may suggest that at temperatures 
not very far above Td, that even at weak coupling the the- 
ory is non-perturbative in the sense that Aq/T is large. 
If so, it is useful then to construct an effective theory in 
terms of the Wilson line L rather than Aq [l|, 0, H2] ■ This 
Lagrangian can also incorporate the global Z(N) symme- 
try for the Polyakov loop. 

As shown in ref. [lj, the electric field in the three- 
dimensional theory for arbitrary Aq is given by 



1 That is, all fields in JTJ are functions of x only, and the action 
is given by an integral of the Lagrangian over space, divided by 
temperature. 



E l {x)^-V(x)D l (x)L(x) 



(4) 



The classical Lagrangian in three dimensions then be- 



comes 



II. THE LATTICE ACTION 
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err 



2^Gl 
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trlLtALl 2 



(5) 



Contrary to sigma models with left-right symmetry, for 
loops there is also a potential. It can be written as an 
infinite sum over all Z(N) neutral loops |13J|. For the 
present purposes, however, we rather write it in terms of 
powers of the fundamental Wilson line. To one loop and 
for constant L 1141. 



L\ 



eff 



-loop 



^E 



1 



ItrL 



n\2 



(6) 



This potential is evidently minimized by the perturbative 
vacuum (L) = 11 (times a phase), for any T. To gener- 
ate eigenvalue repulsion and a phase transition in infinite 
volume, Pisarski suggested [l| to add non-perturbative 
contributions such as 



C?„ =B f T 2 \tYh\ 2 



off 
■'non- pert. 



(7) 



with Bj a "fuzzy" bag constant (see, also, refs. [la fig). 
At sufficiently low temperature, (jTj> dominates over the 
perturbative potential ((6]) and induces a transition to 
a confined phase with (trL) = 0. It was further sug- 
gested in [l| that terms such as (J7]) lead to "repulsion" of 
the eigenvalues of the Wilson line in some temperature 
range above T^. If so, then the distribution of eigenval- 
ues should deviate from a sharp peak near 1 for non- 
asymptotic temperatures. 

In this paper, we perform Monte-Carlo simulations of 
the effective theory given by (0[7]) on a three-dimensional 
lattice. The theory is defined with a spatial cutoff 
on the order of the inverse temperature as ([5]) is non- 
renormalizable in three dimensions and is valid only over 
distance scales larger than 1/T 2 . We shall focus in 
particular on measuring the eigenvalue distribution both 
above and at the (de-)confining phase transition, thereby 
testing the presence of eigenvalue repulsion in the phase 
transition region. We presently employ several approxi- 
mations which simplify the simulations drastically. Most 
importantly, the present simulations neglect the mag- 
netic sector, 



A, = 



(8) 



Hence, the gauge theory is essentially reduced to a sigma 
model. A precise matching of the couplings in the ef- 
fective theory to correlation functions measured in the 
continuum limit of the original four-dimensional theory 
is beyond the scope of this paper. Furthermore, we ne- 
glect all but the n = 1 term in (|6|), which can then be 
combined with the non-perturbative potential IJ7]). 



2 A related renormalizable theory has been formulated in ref. [SJ. 
Ref. [l7j derived the relations between lattice and continuum 
theories to leading order in lattice perturbation theory. 



Our general three-dimensional lattice action includes 
kinetic (nearest-neighbor interaction) and mass terms, 

S = -~/3^tr (i^LJ+h.c.) -m 2 ^|trL,| 2 , (9) 

(i'j) > 

where L denotes SU(2) Wilson lines in the fundamen- 
tal representation, i labels sites, and (ij) labels links. 
We employ periodic boundary conditions. The kinetic 
term is invariant under global SUl(2) x SUr(2) trans- 
formations while the mass term breaks it to SU(2). The 
weak-coupling limit of the original four-dimensional the- 
ory corresponds to large (3. The partition function in- 
volves an integral over the invariant SU(2) measure [dL] 
at each site, 



Z = 



n ^ 



,-s 



(10) 



III. MEAN-FIELD APPROXIMATION 

The mean-field approximation has been discussed in 
detail in refs. [1, [13, Gil 03 • We briefly review the main 
steps and results as required for our present purposes. 

Replace all 2d nearest neighbors of any given site in ([9]) 
by a fixed matrix L, where d = 3 is the number of spatial 
dimensions. This defines a single-site free energy, 



-N*F aa (L) 



,JV 



fill 



where N s is the number of sites per spatial dimension 
and 



Z ss = I [dL] exp d(i tr ( LL f + h.c. 
Consistency requires that 



m 2 ItrLI 2 



<(l; 



w I 



be equal to 



—rr, — = — log Z ss (L) 







(L)? fc = -=- (L)r fc (L) ifc 

o[L)i k 



(12) 
(13) 

(14) 



It follows that L minimizes a mean-field free energy de- 
fined as 







d_ 



^ m r(L) 



■Fmf(L-) = ^ s (L)+d/3trL f L 



(15) 
(16) 



To proceed, we assume that L is proportional to the unit 
matrix, L = Zll (for two colors, £ can be chosen to be 
real), so that 



-^•..(Z) _ 



= /[dL] exp^d/^trQL + L^ + m^trLI 2 ] 



(17) 



The action is a function only of the trace of the integra- 
tion variable, so that we can write 



exp diag (i(f>, —i<p + 2-Kin) 



(18) 



with n an arbitrary integer, and employ Weyl's parame- 
terization 



[dL] — d0 |A(0)| 2 = d<^sin 2 



(19) 



where A(0) denotes the Vandermonde determinant. Up 
to an overall constant then. 



-^•® = / dcos 



exp 



Ad(3£ cos (j) + 4m 2 cos 2 (f> + - log(l 



(20) 



This integral could now be evaluated analytically in a 
saddle-point approximation. However, we have found 
that for d = 3 the analytical result is too inaccurate to be 
useful in practice, in particular in the interesting region 
of j3 and m 2 . Therefore, we have rather tabulated ([20]) 
as a function of I. The expectation value £q of trL/2 is 
then given by the location of the minimum of 



Fmi(t)=f..W + 2d0e' 



IV. RESULTS 



(21) 



A. The model with global SUl(2) x SUr(2) 

symmetry 

We begin with the pure nearest-neighbor interac- 
tion model with no loop potential, corresponding to 



eqs. (|9TT0l) with 



S 



0[18|: 



-£/»£*( 

<«i> 



L,L] 



h.c. 



(22) 



Note that in (f2"2"|) the basic degrees of freedom are the 
Wilson line matrices, or their eigenvalues; the model 
therefore differs from others which deal only with the 
trace of L, such as S ~ — (3J2 (trLj trLt + c.c.) [2Cj |. 
Alternatively, one may consider nearest-neighbor inter- 
actions between Polyakov loops in arbitrary representa- 
tions nm m. 

We expect that for small (3 there is a phase where the 
adjoint fields 



£ a (x) = — trL(x)r a 



(23) 



as well as the singlet field (which is actually the Polyakov 
loop) 



l(x) = -trL(x) 



(24) 
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FIG. 1: The integrated autocorrelation time as a function of 
the coupling /3 for various lattices. 



are massive 3 . Furthermore, the expectation value of the 
"length" of L, 



u = \/trL L/2 , uq = (u) 



(25) 



should vanish also. The bar stands for the average over 
the volume for any given configuration: 



1 v- 



iV 3 



(26) 



while (•) is the average over configurations. Note that 

trL f L/2 = 2(trL/2) 2 - trL 2 /2 =t - \, where £ 2 is 
the Polyakov loop with Z(N) charge two [22| (which is 
neutral when N = 2). 

For sufficiently large j3, on the other hand, the Wilson 
lines at different sites have to align in order to minimize 
the action (f2"2"|) . Hence, for a given configuration (resp. 
Metropolis time) L should be non-zero. However, its di- 
rection in group space will rotate from configuration to 
configuration, implying (L) = 0. To monitor the transi- 
tion to an ordered phase at large j3 we therefore use (u) 
rather than (tr L) as order parameter [181 ] . Alternatively, 
one could add a weak background field, — /itrL, which 
is then taken to zero after the extrapolation to infinite 
volume has been performed. 

The regimes where (u) — and (u) ^ 0, respectively, 
are separated by a second-order phase transition at some 
critical j3 c [18| which we determine numerically. This 
transition is associated with spontaneous breaking of the 
SU L (2)xSU fl (2) symmetry to SUy(2), where three Gold- 
stone modes appear. 

The ensemble average denoted by (•} should be per- 
formed over statistically independent configurations. It 



We assume hermitian generators normalized according to 



tr-r a 



: 28 a 



is therefore necessary to determine the autocorrelation 
time of the Monte-Carlo algorithm as a function of j3. 
This is done via the "rebinning method" [23| as follows. 
First, we group the sequence Oi of measurements of a 



given operator into N bs bins of size Nb, 




a+i)jV(,-i 
o 3 (N b ) = — Yl °*> ( 2? ) 

b i=jN b 


o 
3 


where j = • • • N bs — 1 labels the bins. Hence, Oj(N b ) 





is simply the mean over the measurements belonging to 
the bin j. We then determine the variance of the new 
sequence Oj{N b ): 



2 



1 



N, 



hs 



N bs -1 

E 

3=0 



(0 3 (N b )-(0)Y 



(28) 



(O) denotes the average of O over all configurations. The 
integrated autocorrelation time corresponding to the bin- 
size N b is given by 



t(N b 



'N b 



(29) 



where er 2 denotes the variance of the original sequence of 
measurements. We then plot Ti nt (N b ) versus N b , which 
eventually approaches a flat plateau (up to rapid oscil- 
lations) . This defines r; n t , which is shown in Fig. [1] as 
a function of (3. In the vicinity of the critical point, 
the Metropolis update algorithm displays the well-known 
critical slowing down phenomenon; r; n t diverges in the 
infinite- volume limit. The measurements obtained on 
N a = 24, 36, 48 lattices can be fitted with the form 
Ti n t ~ N a , with the scaling exponent 



v T = 0.72(4) 



(30) 



Away from (3 C the autocorrelation time decreases. Notice, 
however, that it increases with the volume even above the 
critical point, while it exhibits the standard behavior for 

In practice, our simulations are performed as follows. 
The initial configuration of SU(2) matrices is chosen ran- 
domly. We then perform a number of thermalization 
steps which is larger than the autocorrelation time Tint 
(determined beforehand in a pre-run) of the order param- 
eter. Subsequently, measurements are performed in time 
intervals slightly larger than Tint ■ We employ a standard 
Metropolis update [23| where all Nf sites are scanned in 
sequence. Typically, we summed on the order of thou- 
sand configurations for each set of couplings. 

Fig. [2] shows the expectation value of the order pa- 
rameter (|2"5)) as a function of (3 on lattices of various 
sizes. Statistical error bars are smaller than the size of 



4 We take O = u defined in eq. H25I I. 
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FIG. 2: The length of the 0(4)-like order parameter from 
eq. (|25p as a function of the coupling /3 for various lattices. 
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FIG. 3: The derivative of the order parameter from eq. (|25p 
with respect to /3 as a function of the coupling (3 on various 
lattices. 



the symbols. There is, clearly, a order-disorder transition 
at (3 C — 0.9. As expected, finite-size effects are visible 
around the transition point {j3 ~ j3 c ). We have verified 
that uq approaches I for (3 3> 1. 

To estimate the infinite- volume limit of (3 C we proceed 
as follows. We first determine the temperature suscepti- 
bility x{0) = duo/df3, as shown in Fig. [3] The location 
of the maximum defines (3 C for any given lattice size. Ex- 
trapolating linearly to 1/N S — 0, we obtain 



(3 C = 0.942(5) 



(31) 



We have verified that the derivative of the average kinetic 
energy per link E ~ Re(trLJLi + i) with respect to (3 
also peaks at the same value of the coupling, which is 
somewhat larger than the estimate from ref. dSj, who 
employed smaller lattices and lower statistics. 

In Fig. 2] we show two time sequences for the Polyakov 
loop just below and far above (3 C . It is clear that be- 
low the phase transition there are only small fluctuations 
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FIG. 4: Time evolution of the volume-averaged Polyakov loop 
below and far above j3 c \ N s — 12 lattice. 



FIG. 5: The inverse spatial correlation length as a function 
of the coupling j3 for various lattices. 



about 0, which decrease on larger lattices. On the other 
hand, at large (3, the Wilson lines partly align and \l\ is 
far from for long time intervals. However, the above- 
mentioned slow rotation of L in group space (in the ab- 
sence of a background field) is clearly visible. 

Next, we consider two-point matrix-matrix correlation 
functions of the form 



Ci(r) = \±J2 1 ~(trLHr )L(r +r)) 



3 JV| f^ 2 

r.ro 



(32) 



The vector r is allowed to point in any of the three prin- 
cipal directions of the lattice (in the positive direction 
only), over which we average. Also, its length is restricted 
to < N s /2 due to the periodic boundary conditions. 

Having determined the two-point function C(r) and its 
statistical error, we perform a \ 2 fit to the functional 
form 



ClW 



1 



rm^ 



(33) 



to extract the inverse spatial correlation length m^ . The 
fits were restricted to r > 4 (in lattice units) such that 
rVdof ~ 1. 

Fig. \5\ displays m^(/3) for lattices of various sizes. Deep 
in the disordered phase correlations extend only over a 
few lattice sites and m^ is therefore independent of the 
volume. This confirms that both the £ a and the Polyakov 
loop £ are massive. Long-range correlations do develop 
near f3 c and mj drops to nearly zero, up to finite-size 

effects. A fit of the form m^{j3 c ) 

48 lattices only) gives the scaling exponent 



N s 1/l/s (N s = 24, 36, 



v(. = 0.938(5) 



(34) 



Quite clearly, there are massless modes (again, up to 
finite-size effects) even above j3 c and hence toj remains 
small. These observations are in line with the behav- 
ior of the integrated autocorrelation time 7"i n t for (3 > [3 C 



mentioned above. We have also measured the correlation 
lengths for Polyakov loops and for the adjoint t a fields 
via fits of the form (f3"3"|) to the two-point functions 



Ci{r) ~ ^(IW^o + r)) 

f,r 

C- e (r) - £(*(r„)-*(r +r) 

r,r 



(35) 
(36) 



We refrain from showing the results here since they 
closely resemble mf(/3) from Fig. \5\ The fact that the 
correlation length for £ appears to diverge even above j3 c 
is probably due to mixing with the Goldstone modes. 

Finally, we determine the distribution of eigenvalues 
of the Wilson lines. For any given configuration (i.e. 
Metropolis time t), we compute the eigenvalues Ai and 
A2 of the Wilson lines L at each lattice site. We introduce 
their difference and average, 



pi(t,x) 

P2(t,x) 



\Xi(t,x) - X 2 (t,x)\ 
\Xi(t,x) + A 2 (t, a;) I 



(37) 



The ensemble of pi(t, x) defines its probability distri- 
bution Pi (pi), and similarly for P2{p2)- These can be 
turned into effective potentials for the sum and differ- 
ence of eigenvalues, respectively, via 

^eff(pi) = -logPi(pi) , V e s(p 2 )=- log P 2 (p 2 ). (38) 

Figures [6] and [7] depict the probability distributions 
P(pi), P(p 2 ) and the corresponding effective potentials 
for j3 = 1. We have also determined these quantities be- 
low the transition (/3 = 0.5) but obtained very similar 
curves. The potential shows evidence for a logarithmic 
divergence at p x — > (or p 2 — > 1); this is expected as 
the group integration measure leads to logarithmic repul- 
sion of the eigenvalues, compare to eqs. (|19l20p . Aside 
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FIG. 6: The probability distributions of the difference (pi) 
and sum (p%) of eigenvalues of the Wilson line for /3 = 1 
obtained on a N s = 48 lattice. 



FIG. 8: The integrated autocorrelation time at (5 = 1 as a 
function of the coupling m 2 for various lattices. 
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FIG. 7: The effective potentials for the difference (pi) and 
sum (P2) of eigenvalues of the Wilson line for /3 — 1 obtained 
on a N s = 48 lattice. The pure integration measure in terms 
of p2, which is given by log(l — pi) , is shown by the points. 



from the effects of the Vandermonde determinant, how- 
ever, the eigenvalue distribution (or the potential) for 
p-2 = (l/2)|trL| is entirely flat. This is illustrated in 
Fig. [7] which compares the pure Vandermonde potential 
log(l — p!) -1 ^ 2 to the actually measured V(p2)- The flat 
eigenvalue distribution is consistent with the free global 
rotations of L observed above. 



B. Action with SU(2) symmetry 



In this section, we add a mass term for the Polyakov 
loop^ = trL/2 Q, 

S = -^/3j2 tT (l^LJ+h.c.) -m 2 ^|trL,| 2 , (39) 

(ij) « 



which explicitly breaks SU L (2) x SU fl (2) to SU(2), L -> 
fjtLfJ, and also respects the Z(2) symmetry for the 
Polyakov loop, t — ► —I. We study the phase structure as 
a function of m 2 at fixed j3. The order parameter for the 
deconfining phase transition is given by the ensemble and 
volume averaged Polyakov loop (|?|) = (|trL/2|), where 
L is defined in eq. (|26|) 5 . In this section, (|Z|) will also be 
denoted as £q. 

In Fig. [5] we show the integrated autocorrelation time 
for the Polyakov loop in the model (|39|) at /3 = 1, as a 
function of m 2 . It indicates that the transition occurs in 
the vicinity of m 2 ~ 0, where r m t grows with the lattice 
volume (critical slowing down). However, contrary to 
Fig. [T] above the transition point Tint is independent of 
the volume. This confirms to our expectation that long- 
range correlations in Metropolis-time should not appear 
for m 2 ^ 0. A fit of the form r int - n] /Ut to the N s = 24, 
36, 48 data gives 



v T = 1.3(4) , 



(40) 



at to 2 = 0. As before, all subsequent measurements were 
performed with configurations that were separated by a 
time interval of Tint (at least). 

Figs. [H and [TU] show the expectation value of the 
Polyakov loop, and its derivative with respect to the cou- 
pling, in a narrow window about the deconfining phase 
transition. Within errors, we find that the critical cou- 
pling is 



m 2 c = 0.000(2) 



(41) 



The transition in terms of m 2 is evidently rather sharp. 
Nevertheless, the scaling of Tint with the lattice size men- 
tioned above suggests a second-order phase transition in 



Taking the absolute value of t before performing the ensemble 
average is required due to the Z(2) symmetry. 
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FIG. 9: The expectation value of the Polyakov loop as a func- 
tion of the coupling m 2 (at (3 — 1) for various lattices. 



s 



0.1 



#>^- 



ifM 



*********** 



*** * §*+** 



5 :i-:ii;ii 

■ • 







N=48 



-0.12 -0.06 



0.06 0.12 



m 



FIG. 11: The inverse spatial correlation length as a function 
of the coupling m 2 (at (3 — 1) for various lattices. 



X 



10 z 



10' 



10 U 



10" 



10" 



x* 
x i - x r 






N s =12 ^^ 
i N g =48— ♦- 

-0.04 -0.02 0.02 0.04 

2 
m 



FIG. 10: The derivative of the Polyakov loop with respect to 
m 2 as a function of the coupling m 2 (at (3 — 1) on various 



lattices. 



infinite volume, which is consistent with the global Z(2) 
symmetry for the Polyakov loop. This is confirmed also 
by the behavior of the inverse correlation length m^(m ) 
shown in Fig. 111! m^ has been determined by the same 
procedure outlined in eqs. (|32I33[) from the previous sec- 
tion, and appears to vanish at m 2 = 0, N s -too; fitting 
m 5 ~ Ns l/n (to the N s = 24, 
scaling exponent 



36, 48 data) gives the 



^ = 2.28(8) 



(42) 



In the deconfined phase at m 2 > 0, the correlation length 
decreases rapidly to about one (in lattice units). It de- 
creases also as one goes to negative values of m 2 , into the 
confined phase, but less rapidly. There, a weak volume 
dependence remains even from N s = 36 to N s = 48. 

The expectation value of the Polyakov loop is shown 
again in Fig. [T^] over a broader range of m 2 . We also 
compare to the mean-field prediction (only for /3 = 1) dis- 
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FIG. 12: The expectation value of the Polyakov loop as a 
function of the coupling rn 2 at various f3. Monte-Carlo results 
obtained on a N a = 12 lattice are indicated by the symbols. 
The line indicates the mean- field prediction for (3 = 1, shifted 
horizontally by m^ ni — m 2 — 0.94. 



cussed in section UTTl which has been shifted to the right 
by Am 2 = 0.94 to match the data far above the transi- 
tion. Such a shift is expected by analogy to the tadpole 
contribution in a scalar theory, for example. Not surpris- 
ingly, mean-field theory works well for large \m 2 \ <^ 0.5 
(far from the transition, to both sides), when the effective 
masses are large and fluctuations are suppressed. Close 
to the phase transition, critical fluctuations invalidate the 
mean-field approximation. 

The transition becomes extremely sharp when [3 is 
large, switching almost instantly from the confined phase 
to a perturbative deconfined phase with £q ~ 1. This be- 
havior is in line with the discussion in sections U and [Til 
positive m 2 and large (3 corresponds to the weak-coupling 
limit of the four-dimensional theory in the deconfined 
phase. It can also be readily understood from the ex- 
pression (|3"9")1 for the action: at large j3 the Wilson 
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FIG. 13: Eigenvalue distributions for m 2 = 0.15 and /3 = 1; 
N 3 = 48 lattice. 



FIG. 14: Eigenvalue distributions for m 2 = —0.05 and (3 — 3; 
iV s = 48 lattice. 



lines at neighboring sites are forced to align such that 
trL<l4/2 w 1. The potential only determines the direc- 
tion of alignment: when m 2 > 0, the preferred direction 
is the unit matrix (eigenvalue attraction); when m 2 < 0, 
the Wilson lines instead live in the subspace spanned by 
the Pauli matrices (eigenvalue repulsion). 

Fig-EHalso shows that the deconfining phase transition 
is shifted to m 2 > when (3 < (3 C . In this limit the align- 
ment of the Wilson lines is enforced by the upside-down 
potential rather than the nearest-neighbor interaction. 

Fig. [T3] depicts the eigenvalue distribution at j3 = 1 
and m 2 =0.15 which exceeds the critical m 2 for decon- 
fmement (since £q — 0.6) but is still far from asymptotic. 
Here, the perturbative potential ^ is partly cancelled 
by the "fuzzy bag" term ([7]) and the eigenvalue distri- 
butions are rather broad. This result demonstrates that 
the "fuzzy bag" term can generate eigenvalue repulsion 
in the deconfined phase, in the regime ~ C corre- 
sponding to moderately weak coupling in the underlying 
four-dimensional theory. 

In the confined phase at m 2 < and > C the Wilson 
lines fluctuate about the non-trivial vacuum L c = ir^, 
or SU(2) rotations thereof 0, [ija, [2J]; this was shown 
in Fig. [S] already. As expected, the fluctuations dimin- 
ish with increasing 0, see Fig. 1141 They are visible 
mostly in the distribution of the average eigenvalue pi 
while P(p\) is rather sharp. This can be understood 
easily by parameterizing the fluctuations about ir^ as 
L - i diag (e^, -e"^), with « 0. Then, 



Pi = 


<t> 2 

| COS 0| ~ 1- — , 


(43) 


92 = 


sin 0| ~ |0| . 


(44) 



For large m 2 one of course approaches the perturbative 
vacuum, as shown in Fig.[T5j The distribution for p^ = \£\ 
peaks near 1 while that for the difference of eigenvalues is 
broader. In the perturbative regime fluctuations can be 
parameterized as L ~ diag (exp i<f>, exp —i<f>), with w 0. 
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FIG. 15: Eigenvalue distributions for m 2 — 0.8 and = 1; 
iV s = 48 lattice. 



Hence, the fluctuations of p\ = | sin 0| are much bigger 
than those of pi = | cos0| ~ 1 — 2 /2. The fact that 
P{pi) ^ as pi — » 0, and P{p2) — * as p-z — »• 1, is 
again due to the integration measure, see eq. (fH)|) . For 
even larger m 2 , both distributions get sharper and their 
maxima move further towards p\ = and p% = 1, re- 
spectively. In all, far above the transition the eigenvalue 
distributions qualitatively exhibit the behavior appropri- 
ate for the perturbative weak-field regime. 



V. SUMMARY AND CONCLUSIONS 

We have performed Monte-Carlo simulations of an ef- 
fective theory of SU(2) Wilson lines in three dimensions. 
The main purpose of this work was a study of eigenvalue 
repulsion in the deconfined phase. We considered the 
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FIG. 16: Schematic sketch of the phase diagram in the /3 — m 
plane, for infinite volume. The presence of an infinitesimal 
background field — ZitrL is assumed. 



action 



*=4*£ 



tr 



L,L] 



h.c.J - m 2 'S~\ |trLj 



(45) 



(ij) 



without gauge fields, A c [ = 0. The kinetic term exhibits 
a global SUl(2) x SU.r(2) symmetry which is broken 
explicitly to SU(2) by the loop potential. Note that a 
L — » fi^Lf^ transformation changes the eigenvalues of 
L, while L — > Jl^Lf? does not. 

The phase diagram is sketched in Fig. [THJ In the ab- 
sence of a potential, at to 2 = 0, (|45|) is essentially a stan- 
dard spin-model. At small /3 the effective mass of the 
Wilson lines is large and they fluctuate independently 
from site to site. Confinement is realized in a trivial way 
since L — > for each configuration, where L denotes 
the volume-averaged Wilson line. This remains true for 
small | ni |. To dcconfine, a large upside-down poten- 
tial (to 2 > 0) is required to align the Wilson lines to 
the unit matrix. Hence, for small (3 the phase transi- 
tion arises due to the effective loop potential, in a regime 
where SUl(2) x SU,r(2) is broken strongly. 

There is a second-order phase transition at j3 c ~ 0.942 
(and m 2 = 0) where the masses (inverse correlation 
lengths) of the Polyakov loop I = tr L/2 and of the ad- 
joint fields £ a = — itrLr a /2 vanish. This is associated 
with spontaneous breaking of SUl (2) xSUh(2) to SU(2), 
where three Goldstone modes appear. We have confirmed 

that the "length" u 2 = tr L L/2 of L aquires a non-zero 
expectation value for > [3 C . Hence, we expect that a 
weak background field — ft.tr L, h — ► 0, shifts the phase 
boundary to m 2 < as indicated in Fig. [TH 

Very large lattice coupling j3 3> 1 corresponds to 
the extreme weak-coupling limit of the original four- 
dimensional theory; the effective theory can neverthe- 



less confine because it incorporates the global Z(N) sym- 
metry for the Polyakov loop. At large (3 fluctuations 
are suppressed and the Wilson lines are again forced to 
align, this time by the nearest-neighbor interaction (ki- 
netic term). The direction of alignment is determined 
by the loop potential. A standard potential with pos- 
itive curvature (m 2 < 0) is minimized by Wilson lines 
with no singlet component, hence eigenvalues repel and 
the theory confines 6 . On the other hand, an upside-down 
potential (to 2 > 0) leads to h(x) ~ 11 and so to eigen- 
value attraction and deconfinement. For (3 ^$> 1 even a 
weak potential suffices to trigger the locking into (or out 
of) the center of the group. This leads to a sharp transi- 
tion directly to a perturbative deconfined phase without 
eigenvalue repulsion. 

We have measured the distributions of the eigenval- 
ues of the Wilson line in the "weakly-coupled" 7 but non- 
perturbative deconfined phase above (but close to) (3 C . 
They show clearly the emergence of eigenvalue repulsion 
even for "temperatures" (i.e. to 2 ) not extremely close to 
the phase boundary. It is only relatively deep in the de- 
confined phase (to 2 > 1) that the distribution of eigen- 
values peaks near 1, which corresponds to the pertur- 
bative vacuum. These results confirm the suggestion of 
ref. []J that eigenvalue repulsion in the deconfined phase 
does arise at moderately weak coupling (j3 <L (3 C ), due to 
fluctuations of the Wilson lines, provided that the non- 
perturbative "fuzzy-bag" term approximately cancels the 
perturbative loop potential. Such a "fuzzy bag" contri- 
bution in the effective theory makes it possible to reach 
the region of small to 2 in the phase diagram. 

In the confined phase at j3 > (3 C the volume-averaged 
Wilson line L approaches the center-symmetric vac- 
uum 0, eh m 



L c = diag (l,z,z 2 



,N-1\ 



(* 



_ 2-ki/N 



), (46) 



which for two colors corresponds to L c = i T3 (up to 
an overall SU(2) rotation). This is due to the fact that 
the Wilson lines align at large j3, and to 2 < favors 
a direction orthogonal to unity. We repeat that this is 
not the case when is small, where instead L — > for 
to 2 ~0. 
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6 We expect that the phase boundary is shifted from m 2 = to 
some smaller value if an infinitesimal background field is applied. 

7 This terminology is motivated by the fact that for /3 > /3 C the 
alignment of Wilson lines is driven by the kinetic term. By con- 
trast, for j3 <S /3c the deconfining phase transition occurs at large 
m 2 > and is driven by the loop potential. 
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presented here were performed at the Center for Scien- 
tific Computing (CSC) at Frankfurt University. Our code 
is based in part on the MILC collaboration's public lat- 



tice gauge theory code , see 

http://physics.utah.edu/~detar/milc.html 
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